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Preliminary study of Bb parameter using Lattice NRQCD 
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We present the preliminary result of the calculation of Bb parameter using lattice NRQCD. The calculation is 
performed on a 16 3 x 48 quenched lattice at p = 5.9 with clover light quark. The use of lattice NRQCD enables 
us to investigate the 1/rriQ correction from the static limit. The observed mass dependence is well described by 
the vacuum saturation. 
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1. Introduction 

The theoretical determination of the B meson 
decay constant fg and the bag parameter Bb are 
necessary to constrain |V^Vt<j|. The lattice tech- 
nique has been considered to be the most reliable 
approach, and a number of authors have calcu- 
lated these quantities in different approaches. In 
recent years, the value of fs is almost being set- 
tled (but still within the quenched approxima- 
tion) [0. On the other hand, Bb has been calcu- 
lated only under the approximation of mb — » oo 
or a -1 3> rrib, and there has been no systematic 
study of its heavy quark mass dependence. In this 
talk we study the 1/toq correction to the static 
limit using the NRQCD action, which is the first 
step to obtain Bb with the realistic 6-quark mass. 



2. Method 

Since the chiral symmetry is broken on the lat- 
tice, several lattice operators contribute to Bb 
through 

B B {m b )=yZ l Bf\ (1) 
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and operators are given by 

6 l L = 6^(1 - 7s)d 6^(1 - 75 )d 

& N = 2b^{l~ l5 )db^(l+ l5 )d 

+ 46(1 - 75 )d 6(l + 75 )d 
6 l s = b(l-j 5 )db(l- l5 )d. 

The operators other than O l l appear with 0(a s ) 
mixing coefficients. Since the one-loop calcula- 
tion of the coefficients Z\ has not yet been done 
for the NRQCD action, we use those in the static 
approximation instead in this preliminary study. 
Z[s at /3=5.9 are given in Table @. 

The method of the calculation of B\ at, s is stan- 
dard. We place a point source of the heavy and 
light quarks at the origin, where the 4-fermi 
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Table 1 

Matching factors at /3=5.9. 
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Size 


16 3 x 48 (t = [-24, 23]) 


Gauge 


Wilson plaquette action 

quenched 100 config. 

/3=5.9, fixed to Coulomb gauge 


Light quark 


SW-Clover action(c sw = 1/uq) 
local source at t = 
KLM norm.|] 


Heavy quark 


Tadpole improved NRQCD 
0(l/m) and 0(l/m 2 )§ 
local source at t = 



Table 2 

Actions and simulation parameters 

operators are constructed, and B and B mesons 
propagate in the opposite direction to each other. 
The ground state extraction is rather easier for 
NRQCD than for the static approximation, even 
though we use the local sink. This is another 
advantage of introducing the 1/toq corrections. 

The actions and lattice parameters are summa- 
rized in Table |. The inverse lattice spacing is de- 
termined from the rho meson mass as a -1 = 1.62 
GeV. For the heavy quark, we performed two sets 
of calculations with X/rriQ and 1/toq NRQCD ac- 
tions as in Ref. |l , in order to estimate the higher 
order effects in the 1/toq expansion. 

3. Results 

In Figure |l| we plot the 1/Mp dependence of 
B\ at, s, with Mp the pseudoscalar heavy-light me- 
son mass. Since B l £* should agree with B 1 ^ with 
infinite statistics, we take an average of B 1 ^ and 
Bp^" in our simulation. Bjy and Bg have signifi- 
cant slopes in 1/Mp, while B^p are almost flat. 
The curvature in 1/Mp, on the other hand, is 
small in any case, and the difference between the 
two different NRQCD actions is negligible. An 
extrapolation to the static limit does not seem 
to agree with the result of || (at /3=6.0) for Bn 
and for B$, although it is not quite clear whether 
this is an inconsistency, with our large statistical 
errors. 

This 1/Mp dependence can be roughly under- 
stood by the factorization hypothesis (or vacuum 
saturation) as follows. If the factorization is exact 

(fad 

on the lattice, B L R becomes unity by definition 
(solid line in Figure |I|). For B^g, after a little 
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Figure 1. The heavy quark mass dependence of 
B' at . The circles are for 0(l/m,Q) calculation and 
the triangles for 0(l/rag). Opaque symbols are 
obtained by interpolations to the physical B me- 
son mass, respectively. Diamonds are the results 
obtained in Ref.0. 
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algebra, we obtain 
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where 5fp^ is the 0(l/mg) correction to the de- 
cay constant fp. 5f B is known to be rather 
large, and as an estimate we use our previous 
calculation with NRQCD Substituting the 
value of Sfg/fs to Eq.(|^) and (||), we obtain 

open squares in Fig.[l] for B^g . The agreement 
between the full calculation and the factorization 
is remarkable. Not only their value itself, but also 
the slope in 1/ toq seem to be well described with 
the factorization hypothesis. 

The 1/Mp dependence of B B (mb) is shown in 
Figure^, which shows a little negative slope. This 
is because of a large cancellation of 1 /Mp depen- 
dence between those in B^ and B$- In the lead- 
ing log approximation, with Kq CD =Q.2 GeV, we 
obtain the scale invariant Bp parameter 

B Bd = (a s {m b ))~ 6 / 23 B B (m b ) = 1.17(9)(4)(3). 

The errors are statistical, higher order in pertur- 
bation theory, and the higher order in 1/toq ex- 
pansion, respectively. For other quantities, we ob- 
tain B Bs /B Bd = 1.00(2) and f B Jf Bd = 1-16(3). 
Some additional systematic uncertainties, such as 
0(a 2 ), 0(a s /m), 0{aa s ) and quenching errors 
still remain to be estimated. 

4. Summary 



To summarize, we find a large 1/toq effects 
in the matrix elements B 1 ^ and B l g l , although 
the 0(1/ nig) correction seems reasonably small. 
Their large mass dependence can be roughly un- 
derstood by the factorization hypothesis (or vac- 
uum saturation). In the present analysis, we 
combine lattice simulation for finite heavy quark 
mass with matching coefficients in the static limit, 
which is still unsatisfactory. For the complete un- 
derstanding of I/ttiq dependence, matching coef- 
ficients with finite heavy quark mass are abso- 
lutely necessary. 
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Figure 2. The heavy quark mass dependence of 
Bp(mb). A diamond is the result of reanalysis || . 
For other symbols, see the previous figure cap- 
tion. 
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